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A Local One-Zone Model of MHD Turbulence in Dwarf Nova Disks 

Takayoshi Sano 1,2 and James M. Stone 1 ' 2 

ABSTRACT 

The evolution of the magnetorotational instability (MRI) during the transition from 
outburst to quiescence in a dwarf nova disk is investigated using three-dimensional 
MHD simulations. The shearing box approximation is adopted for the analysis, so 
that the efficiency of angular momentum transport is studied in a small local patch 
of the disk: this is usually referred as to a one-zone model. To take account of 
the low ionization fraction of the disk, the induction equation includes both ohmic 
dissipation and the Hall effect. We induce a transition from outburst to quiescence 
by an instantaneous decrease of the temperature. The evolution of the MRI during 
the transition is found to be very sensitive to the temperature of the quiescent disk. 
As long as the temperature is higher than a critical value of about 2000 K, MHD 
turbulence and angular momentum transport is sustained by the MRI. However, MHD 
turbulence dies away within an orbital time if the temperature falls below this critical 
value. In this case, the stress drops off by more than 2 orders of magnitude, and 
is dominated by the Reynolds stress associated with the remnant motions from the 
outburst. The critical temperature depends slightly on the distance from the central 
star and the local density of the disk. 



Subject headings: accretion, accretion disks — diffusion — instabilities — MHD 
novae, cataclysmic variables — turbulence 



1. INTRODUCTION 

Dwarf novae are close binary systems that consist of a white dwarf and a Roche-lobe-filling 
secondary star. The matter overflowing from the secondary star forms an accretion disk around 
the white dwarf. Dwarf nova systems show repetitive outbursts in which the luminosity of the 
disks increase rapidly. The duration of the outbursts ranges from a few to 20 days, and the 
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recurrence time is about 20 - 300 days. The disk instability model (Osaki 1974), in which the 
evolution of the disk is regulated by changes in the rate of angular momentum transport, is 
generally accepted as the explanation of the outbursts (Cannizzo 1993; and references therein). 
Since ordinary molecular viscosity is too low to explain the evolutionary timescale of the disks, an 
anomalous stress associated with turbulence and characterized by the a parameter (Shakura h 
Sunyaev 1973), is thought to be the source of angular momentum transport. Detailed comparisons 
between observed light curves of dwarf novae and the theoretical disk models suggest that the 
viscous parameter a must vary between the hot (outburst) and cold (quiescent) state, and the 
amplitude of a^ot an d a coid is typically of the order of 0.1 and 0.01, respectively. However, while 
the origin of the anomalous stress has become better understood in recent years (see below), the 
reasons for the difference between a^ot an d a C oid are still unclear. 

The magnetorotational instability (MRI; Balbus &; Hawley 1991) is the most promising source 
of the anomalous stress. MHD turbulence driven by the MRI can transport angular momentum by 
the Maxwell (magnetic) stress. During the hot state, the disk gas is fully ionized so that the ideal 
MHD approximation is appropriate. Numerical simulations that adopt ideal MHD have shown 
that the Maxwell stress caused by the MRI can give a ~ 0.1 to 0.01, where the saturation level is 
determined by the geometry and strength of the magnetic field (e.g., Hawley, Gammie, & Balbus 
1995; 1996). Therefore, MHD turbulence can be the primary mechanism of angular momentum 
transport at least during outbursts, that is it can account for cchot- 

During the cold state, on the other hand, the gas is only weakly ionized. Gammie & Menou 
(1998) have pointed out that ohmic dissipation could modify the nature of the MHD turbulence 
at quiescence, and this could make the difference between ahot and a co id- Calculations of the 
ionization fraction in dwarf nova disks reveal that the Hall effect as well as ohmic dissipation must 
be considered if the temperature is T <J 2000 K (Sano Sz Stone 2002a). At low temperatures, 
ohmic dissipation can suppress the growth of the MRI (Jin 1996; Sano Sz Miyama 1999). If 
turbulence dies away completely as a result of ohmic dissipation, the Maxwell stress cannot 
provide the required anomalous stress during quiescence (Menou 2000). Thus, it is important to 
calculate numerically the amplitude of the Maxwell stress at the cold state. 

In this paper, we examine the behavior of MHD turbulence during the transition from an 
outburst to quiescence using local three-dimensional MHD simulations. These simulations include 
the dominant non-ideal MHD effects as determined from a self-consistent calculation of the 
ionization state in the disk. Since we adopt a local approximation to study the turbulent stresses 
(usually referred to as a one-zone model, e.g., Mineshige Sz Osaki 1983; Cannizzo Sz Wheeler 
1984), our calculations do not include global effects (such as spiral shocks in the disk, Sawada, 
Matsuda, Sz Hachisu 1986) which may contribute to angular momentum transport. 

The plan of this paper is as follows. In § 2, we calculate the ionization fraction in dwarf 
nova disks by solving the Saha equation to estimate the magnitude of nonideal MHD effects at 
quiescence. Our numerical method is described in § 3. The results of numerical simulations which 
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examine the dependence of the stress on the geometry of magnetic field and the decay timescale of 
turbulence are presented in § 4. In § 5, we discuss the activity of MHD turbulence and the source 
of angular momentum transport during quiescence. 



2. IONIZATION STATE OF THE DISKS 

During outburst, the temperature of the disk is above 10 4 K, so that the gas is fully ionized 
and the ideal MHD approximation holds. In this case, the disk will be dominated by MHD 
turbulence driven by the MRI, with the plasma beta j3 ~ 10 2 - 10 4 (e.g., Hawley, Gammie, & 
Balbus 1995; 1996), so that the Alfven speed is typically v& ~ 10 4 - 10 5 cm s — 1 . During the 
transition to quiescence, the temperature drops to a few thousands K. The ionization fraction 
decreases dramatically, and thus nonideal MHD effects must be included. The induction equation 
is then 

FIR r AirnJ T y R ( T y Rl y Rl 

(1) 
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where v is the neutral velocity, 77 is the magnetic diffusivity, J is the current density, c is the speed 
of light, e is the electron charge, n e is the number density of electrons, 7, is the drag coefficient, 
and pi and p are the mass densities of ions and neutrals, respectively. Each term in the brackets 
of equation (|T|) from left to right means the inductive term (I), ohmic dissipation (O), Hall effect 
(H), and ambipolar diffusion (A). 

The importance of the nonideal effects can be measured by the relative ratio of each term in 
equation (1) to the inductive term. These ratios are determined mainly by the ionization fraction. 
At a typical density of dwarf nova disks (n n ~ 10 18 cm -3 , where n n is the number density of 
neutrals), ohmic dissipation and Hall effect are very important and affect the evolution of the 
MRI (Sano & Stone 2002a). The size of the ohmic dissipation is given by the magnetic Reynolds 
number; 

Reu = I , (2) 

where f2 = (GM/r 3 )" 1 / 2 is the angular velocity, G is the gravitational constant, M is the mass of 
the central object, and r is the distance from the central star. The magnetic diffusivity is given by 

, = ^ = 234 f V/ 2 cmV\ (3) 
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where v e = n n (av) e is an effective collision frequency for electrons. Here we take 
(av) e = 8.28 x lO^^T 1 / 2 (Draine, Roberge, k Dalgarno 1983). The magnetic Reynolds 
number ReM is the ratio of the inductive term to the ohmic dissipation [Reu = I/O; Sano &; 
Stone 2002a), and a smaller ReM represents a more resistive case. The Hall effect on the MRI is 
characterized by the Hall parameter (Balbus & Ter quern 2001); 

cBQ 



X = t , (4) 
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which is proportional to the ratio between the Hall and the inductive term (X = 2H/I; Sano & 
Stone 2002a). If the Hall parameter X is large, the Hall term is important. 

Both the magnetic Reynolds number and Hall parameter can be written as a function of the 
number density of electrons and the temperature; 



ReM = 68 ' s 
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- W J , (5) 
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Interestingly, the product ReMX is independent of the ionization fraction; 

The number density of electrons is a function of the temperature, because the free electrons 
in dwarf nova disks are mainly generated by thermal ionization of Al and K. Figure |] shows the 
temperature dependence of the electron abundance n e /n n . We have solved the Saha equation 
assuming Solar abundance for the case of n n = 10 18 cm -3 , which is a typical density of dwarf 
nova disks. The ionization fraction decreases as the temperature decreases. As seen from the 
figure, the fraction is very sensitive to the temperature when T <J 2000 K. For example, the 
abundance is n e /n n « 3.5 x 10~ 8 at T = 2000 K, but it is reduced by 2 orders of magnitude 
(n e /n n « 3.7 x lO" 10 ) at T = 1500 K. 

Figure || also shows the temperature dependence of the parameters ReM and X, which are 
given by equations (||) and (||), respectively. Here we assume va = 10 5 cm s _1 , M = M & , and 
r = 10 10 cm. When ReM ^ 1 or X > 1, these nonideal MHD terms can modify the linear growth 
and nonlinear evolution of the MRIH (Sano, Inutsuka, & Miyama 1998; Wardle 1999; Fleming, 
Stone, & Hawley 2000; Balbus & Terquem 2000; Sano & Stone 2002a,b). Then, nonideal MHD 
effects are clearly inefficient if the temperature is higher than 3000 K (ReM S> 1 and X >C 1), 
so that ideal MHD may be valid. However, both the ohmic dissipation and Hall effect become 
important when T 2000 K (ReM ^ 1 and X ^ 1). Therefore, these nonideal MHD effects 
must be included in the analysis when the temperature of the disk becomes below 2000 K. At 
this range of the temperature, the electron abundance and number density is n e /n n <^ 10 -8 and 
n e <^ 10 12 cm -3 , respectively. 



3 The definition of the magnetic Reynoids number given by equation (g) is different from that used in the previous 
work (Re' M = c s H/ri; Gammie & Menou 1998; Menou 2000), and thus the critical value is also different (see Sano & 
Stone 2002a,b). 
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For a typical condition of dwarf nova disks, the relative ratio of the ambipolar diffusion to the 
inductive term is given by 

a n _ 6 / ni y 1 ( m \ 1/2 ( r Y 3/2 

T"^" 6 ' 1Xl ° Uo 12 cm^J [m^J vlOi^J ' (8) 

where we assume that 7, = {av)i/(rrii + m n ), (av)i = 1.9 x 1CP 9 cm 3 (Draine et al. 1983), the ion 
and neutral particle mass are rrii = 30mH and m n = 1.27m-H, and mn is the mass of a hydrogen 
atom. As mentioned in Gammie and Menou (1998), this ratio given by equation (||) is so small 
that the ambipolar diffusion does not affect the evolution of the MRI (Blaes & Balbus 1994; 
Brandenburg et al. 1995; Hawley & Stone 1998). Therefore, we ignore the ambipolar diffusion in 
the numerical simulations in this paper. 

3. NUMERICAL METHOD 

We study the evolution of the MRI and MHD turbulence during the transition to quiescence 
using three-dimensional MHD simulations. The shearing box approximation (Hawley et al. 1995) 
is adopted to investigate the local evolution of the instability. The equations of adiabatic MHD 
with the ratio of specific heats 7 = 5/3 are solved in a Cartesian frame of reference (x, y, z) 
corotating with the disk at the angular velocity f2, where x is the radial, y is the azimuthal, and 
z is the vertical direction. The induction equation includes terms that describe ohmic dissipation 
and the Hall effect. The details of the equations and our numerical scheme are described in Sano 
& Stone (2002b). Vertical gravity is ignored, so that the calculations begin with a uniform density, 
p = po, and a uniform pressure, P = P$. The magnetic field structure of dwarf nova disks is 
still highly uncertain. If the white dwarf has a strong magnetic field, then the dipole field may 
penetrate the disk vertically. On the other hand, a toroidal field amplified by the differential 
rotation may be dominant in the disk. In this case, the net flux of the vertical field could be 
negligible. Nonlinear evolution of the MRI started with a toroidal field is quite similar to that 
with a zero net flux B z (Sano & Stone 2002b). Therefore, in this paper, we consider two kinds 
of initial field geometries; a uniform vertical field, B z = Bq, and a zero net flux vertical field, 
Bz = Bq sm(27rx/H), where H = (2/^) 1 I 2 c s q/Q, is the scale height of the disk and c s q is the initial 
sound speed. We choose normalizations with po = 1, H = 1, and Q = 10~ 3 . The magnetic field 
strength is characterized by the plasma beta /?o = Po/(Bq/8it) = 3200. This is about 100 G in 
dwarf nova disks, which corresponds to a dipole field of a white dwarf with a moderately strong 
magnetic field (~ 10 6 G) at a radius r = 10 10 cm. The size of the simulation box is H x AH x H 
for all the models calculated in this paper. The initial ratio of the MRI wavelength, 2itvao/Q, to 
the box height is about 0.11. A standard grid resolution of 32 x 128 x 32 uniform zones is used for 
most of our calculations. 

To set up a disk in outburst, the induction equation is integrated for 50 orbits assuming ideal 
MHD. The simulations are then restarted with nonideal MHD effects included in the induction 
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equation, and the models are evolved for a further 50 orbits. Although the transition should 
proceed at the thermal timescale of the disk, our technique assumes an instantaneous change of 
the temperature at 50 orbits, which is thereafter kept constant. The quiescent disk is characterized 
by the non-dimensional parameters Rcm and X. We have computed 4 different models with 
(ReM,X) = (100,0.01), (10,0.1), (1,1), and (0.1,10), which correspond to the temperature 
T = 2800, 2200, 1800, and 1600 K, respectively (see Fig. 1). 

As seen from equations (||) and @ , the parameters Re m and X are functions of the magnetic 
field strength. Because the vertical component of the field is essential for the evolution of the MRI 
(Sano &; Stone 2002b), we define these parameters (ReM,50 and A50) using the volume-averaged 
vertical field strength (B 2 ) 1 / 2 and squared Alfven speed {v\ z ) at 50 orbits, i.e., ReM,50 = ( v Az) A?^ 
and X 50 = c(B 2 z ) l l 2 Sl/2-nen e (v\ z ). Then the size of r] and n e after 50 orbits are obtained from 
Rcm,50 and A50, respectively. For simplicity, the electron abundance n e /n n is assumed to be 
constant in time and uniform in space throughout the calculation. Since the pressure has little 
effect on the evolution of the MRI, we find it does not matter whether or not the pressure is 
changed to be consistent with the reduced temperature of the disk after 50 orbits. We have not 
changed the pressure in the models reported here. 



4. NUMERICAL RESULTS 

4.1. Models without a Net Flux of the Magnetic Field 

First, models started with a zero net flux B z are studied in this subsection. All the calculated 
models are listed in Table |l[ To obtain MHD turbulence in the hot state, model S is calculated 
for 50 orbits using the ideal MHD equations. The models of quiescent disk are labeled SA, SB, 
SC, and SD. For each model, ohmic dissipation and the Hall term are included after 50 orbits. 
The size of the parameters Rcm, 50 and X50 are listed in the Table. The magnetic diffusivity and 
the electron abundance at quiescence are calculated using the magnetic field strength of model S 
at 50 orbits. The turbulent stress is the sum of the Maxwell stress, wm = —B x B y /4:ir, and the 
Reynolds stress, wr = pv x 5v y , where 5v y is the perturbed azimuthal velocity from the Keplerian 
flow. The a parameter of Shakura & Sunyaev (1973) is given by a = (wm + wr)/Pq. The time- 
and volume-averaged stresses (((wm)) and {{wr))) are listed in Table [|. The time averages are 
taken from 75 to 100 orbits. 

Figure [| shows the time evolution of the volume-averaged Maxwell stress {wm)/Pqi where 
time is measured in orbits i ro t = 27r/fl. The ideal MHD model (S) is shown by the dashed curve. 
The saturation level of the Maxwell stress in model S is of the order of 0.01 throughout the 
evolution. The solid curves after 50 orbits are for the models that include the nonideal MHD 
effects. From top to bottom, they are models SA (T = 2800 K), SB (2200 K), SC (1800 K), and 
SD (1600 K). 
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The suppression of the MRI is mainly caused by ohmic dissipation (Sano & Stone 2002b). 
When ohmic dissipation is inefficient (.ReM,50 > 100), the amplitude of the MHD turbulence 
maintains a level similar to the ideal MHD case. At ReM,m = 10, the Maxwell stress is reduced 
by a factor of 4, but still has a significant amplitude {{{wm))/Po = 4.7 x 10~ 3 ). However, if 
Rcm,50 ^ 1 or T ^ 1800 K, the MHD turbulence dies away very quickly. As the temperature 
decreases, the amplitude of the magnetic stress decreases dramatically. For the model SD, the 
magnetic field and stress are almost vanished at the end of the calculation, i.e., 50 orbits after the 
transition. The duration of quiescence in dwarf novae is typically more than a few tens of days, 
which corresponds to over 1000 orbits. Thus, the magnetic field in the disks can disappear locally 
during quiescence if the temperature is less than about 1500 K. 

Figure || shows the saturation level of the Maxwell ((%)) and Reynolds stress ((wr)) as a 
function of the magnetic Reynolds number ReM,50- The stresses are normalized in terms of the 
initial pressure Pq, so that they have the same dimension as the a parameter. The upper and 
lower arrows indicate the Maxwell and Reynolds stress for the ideal MHD model {R&m —* oo). As 
shown in Figure ^, active MHD turbulence can be sustain by the MRI if the magnetic Reynolds 
number is much larger than unity. Angular momentum transport is dominated by the Maxwell 
stress, which is about 4 times larger than the Reynolds stress. 

When ReM, 50 ^ 1> the MRI is suppressed by the ohmic dissipation and MHD turbulence dies 
away immediately. The saturation level of the Maxwell stress is quite small (((wm)} / Po ^ 10~ 6 ) 
compared with the less resistive runs (S, SA, and SB). This dependence on the magnetic Reynolds 
number is consistent with the results of previous work (Fleming et al 2000; Sano & Stone 2002b). 
While the Maxwell stress decreases as the temperature decreases, it is found that the Reynolds 
stress has an amplitude which is comparable to or larger than the magnetic stress. Since we do 
not include an explicit kinematic viscosity in our calculations, and since the numerical viscosity 
of our algorithm is much lower than the explicit magnetic diffusivity added to model quiescence, 
epicyclic oscillations remain even after the magnetic field disappears. Although the amplitude of 
the Reynolds stress is very small (a ~ ((wr))/Po ~ 10~ 6 ), these motions, which are remnants of 
the turbulence during outburst, could give non zero accretion stress during quiescence. However, 
such motions may be swamped by global effects which cannot be studied in our local simulations. 

The reduction of the Maxwell stress due to the nonideal MHD effects is independent of the 
numerical resolution. We have performed the same simulations as shown by Figures - I with 
double the resolution (64 x 256 x 64). The saturation level of the stress in the high-resolution 
runs are listed in Table [TJ. For the ideal MHD case, the Maxwell stress in the high-resolution 
run is in fact slightly smaller than that in the standard run. For less resistive models (SA and 
SB), the Maxwell stress is always larger than the Reynolds stress and the saturation level of the 
stress is almost unrelated with the magnetic Reynolds number. When Re m ^ 1, on the other 
hand, the ohmic dissipation suppresses the MRI and MHD turbulence. The Maxwell stress is 
((wm))/Pq = 1-0 x 10 -6 and 1.0 x 10~ 10 for model SC and SD, respectively. For these cases, 
angular momentum transport is dominated by the Reynolds stress with a similar amplitude to the 
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standard-resolution runs (a ~ 10 6 ). 

4.2. Models with a Net Flux of the Vertical Field 

Next we consider the effect of the initial field geometry. Models started with a uniform B z 
are listed in Table 0, along with the model parameters and saturation level of the Maxwell and 
Reynolds stress. With this field geometry, the saturation level of a in the ideal MHD run (Z) is of 
the order of 0.1. This is an order of magnitude larger than that in the zero net flux B z model (S). 
For models with a net flux of vertical field, a two-channel flow forms repeatedly at the nonlinear 
phase of the MRI (Sano & Inutsuka 2001). The channel flow is destroyed through magnetic 
reconnection triggered by the parasitic instability (Goodman & Xu 1994). Since the magnetic 
energy and stress increase exponentially when the channel flow is growing, the time evolution of 
the magnetic stress shows frequent spike-shaped excursions. The growth of the two-channel flow 
can be seen only in the uniform B z runs, and this causes the relatively larger magnetic stress 
compared with the zero net flux B z runs (Sano & Stone 2002b). 

For models of a quiescent disk (ZA, ZB, ZC, and ZD), we choose the same sets of ReM,50 and 
X50 as for the zero net flux B z runs. In Figure |||, the Maxwell and Reynolds stress are depicted by 
the filled and open circles, respectively. The dependence of the stresses on the size of the nonideal 
effects is qualitatively quite similar to that for the zero net flux B z runs. When the magnetic 
Reynolds number is larger than 10, the stress is almost independent of ReM,50 and X^q. The 
difference in the Maxwell stress between the ideal MHD (Z) and the ReM,50 = 10 (ZB) runs is 
only 10 %. The Maxwell stress always dominates the Reynolds stress by a factor of 5. 

If ReM,50 <^ 1, on the other hand, the Maxwell stress becomes inefficient (((wm))/Po ~ 10 -6 ) 
as a result of the suppression of the MRI. For models ZC and ZD, the spatial fluctuations in the 
magnetic field disappear soon after 50 orbits and the field geometry returns to a uniform B z . 
The a parameter is dominated by the Reynolds stress, and the saturation level of {{wr))/Po is 
larger than those in the zero net flux B z models (SC and SD). We also calculated a model with 
Rzm,50 = 0-01 and X^q = 100 (T = 1400 K). In this case, the critical wavelength for the MRI is 
longer than the scale height of the disk (i.e., the simulation box size), because ohmic dissipation 
stabilizes the shorter-wavelength fluctuations. The saturation level of the Maxwell and Reynolds 
stress are ((wm))/Po = —3.6 x 10 -6 and ((wr))/Po = 2.6 x 10 -3 , respectively. The Reynolds stress 
still has the largest amplitude, but this is a remnant of the turbulence during outburst. Once 
again, we emphasize the amplitude of the Reynolds stress may be different in global calculations. 

4.3. Decay Timescale of MHD Turbulence 

When the magnetic Reynolds number is less than unity, MHD turbulence immediately dies 
away after the transition. We define the decay time tdec of the turbulence as the time it takes 
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for the magnetic energy to fall to 1/e of its value at 50 orbits. The decay time for models 
with ReM,50 = 1 (SC and ZC) is t^ cc /t rot = 0.36 and 0.12, where i rot is the orbital time. The 
Rcm,50 = 0.1 runs show much faster decay; tdecArot = 0.028 and 0.0077 for models SD and 
ZD, respectively. The decay time in the high-resolution runs is slightly shorter than that in the 
standard runs (idecArot = 0.16 for model SC and idecArot = 0.012 for model SD). Thus all the 
decay times we obtained are less than an orbit. 

The dissipation timescale of the magnetic field is given by tdis = L 2 /rj, where L is a typical 
length scale. The turbulent motion is driven most effectively at the most unstable wavelength for 
the MRI (A ~ v\/Q), which is much smaller than the disk scale height. We find that the decay 
times idee obtained by the numerical simulations are equivalent to the dissipation timescale with 

l = v A /n-, 

t rot V 27T 

This relation is valid only when the magnetic Reynolds number is less than unity. 



5. DISCUSSION 

5.1. The Critical Temperature and Transition Timescale 

When the magnetic Reynolds number falls to unity, MHD turbulence is suppressed by ohmic 
dissipation independent of the amplitude of the Hall term (Sano & Stone 2002a, b). The critical 
temperature corresponding to Reu = 1 is about 2000 K for our model. As seen from equation (|5|), 
however, the critical temperature depends on the neutral density n n and the distance from the 
central star r. Figure || shows the critical temperature T CI .; t as a function of r for the cases of 
the neutral density n n = 10 17 , 10 18 , and 10 19 cm -3 . We assume the Alfven speed va = 10 5 cm 
s _1 and the mass of the central star M = Mq. At a given temperature, the electron abundance 
decreases as the neutral density increases, so that the critical temperature is slightly higher when 
the density is higher. The MRI wavelength (~ VA_/tt) is shorter at a smaller radius, because 
the angular velocity is larger. Therefore, the MRI can be suppressed with a smaller amount of 
diffusivity r/ in the inner regions of the disk. The critical temperature is about 2500 and 1500 
K at r = 10 9 and 10 11 cm, respectively. Since the magnetic Reynolds number, or the electron 
abundance, has a steep dependence on the temperature, the range of the critical temperature is 
quite narrow and the difference is at most a factor of 2 over the entire disk. 

In our numerical analysis, we assume a fixed temperature during quiescence. Thus the electron 
abundance at quiescence is assumed to be constant. In a real disk, the critical temperature could 
be a little higher due to a runaway decay of MHD turbulence (Menou 2000). The decrease of the 
temperature makes both the electron abundance lower and nonideal MHD effects more important. 
The turbulent stress is therefore suppressed further, so that the temperature falls even more. 
However, a very weak dependence of the stress on the temperature can be seen in our models when 
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T ^ 3000 K. This indicates that the runaway temperature should be less then 3000 K. Although 
a self-consistent study of the evolution of the temperature and electron abundance is needed to 
obtain a real T cr j t , Figure [5] may give a fairly good estimate of the critical temperature in dwarf 
nova disks. 

We also assume the transition to quiescence occurs instantaneously. In actual systems, 
however, this may proceed at the thermal timescale ithj where ith/^rot ~ l/ a <i 100 (Frank, King, 
h Raine 1992). Suppose the temperature decreases gradually. As long as the temperature is 
higher than the critical value, we have shown that the accretion stress is nearly constant. But, 
once the temperature drops below this critical value, the MRI is suppressed, the turbulence decays 
within an orbit, and the a parameter drops by more than 2 orders of magnitude. Therefore, even 
though the decrease of the temperature is gradual, the transition timescale of the stress may be 
very much shorter. 

5.2. Onset of the Next Outburst 

Dwarf nova systems undergo recurrent outbursts. The matter from the secondary star 
accumulates in the disk during quiescence. When the surface density exceeds a critical value at 
some region in the disk, the thermal instability sets in and the temperature goes up. Here we 
examine the behavior of the MRI at the transition to the hot state. 

Figure |6| shows the time evolution of the Maxwell and Reynolds stress for model ZC until 100 
orbits. After 100 orbits, the induction equation for the ideal MHD is used in order to imitate the 
onset of the next outburst. As seen from the figure, the growth of the MRI starts immediately at 
the transition to the outburst and the magnetic field is amplified exponentially. In a few orbits, 
the Maxwell stress reaches the same saturation level as that before 50 orbits. The key to this 
evolution is the existence of a net vertical flux. During quiescence (from 50 to 100 orbits), MHD 
turbulence is suppressed, but the magnetic flux within the shearing box is conserved throughout 
the evolution. Thus the field geometry at 100 orbits is a nearly uniform vertical field. When the 
nonideal MHD effects are turned off, the linear growth of the MRI starts and MHD turbulence is 
initiated. 

This evolution may be different if the system has no net magnetic flux. For example, if the 
magnetic field is completely dissipated during quiescence (e.g., model SD), MHD-driven activity 
cannot occur during the next outburst. However, if enough magnetic field remains for there to be 
unstable modes even in a small region of the disk, MHD turbulence develops in that region and 
finally may spread over a large part of the disk (Hawley &; Balbus 1991). Therefore, it is quite 
important to understand the global structure of the magnetic field during quiescence, including 
the magnetosphere of the white dwarf and the magnetic field supplied by accreting material from 
the secondary star (Meyer & Meyer-Hofmeister 1999). 
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5.3. Beyond the Local Model 

Since we focus on the local behavior of MHD turbulence in this paper, neither global 
structural effects nor global instabilities have been considered. In real disks, global effects such as 
spiral shocks (Sawada et al. 1986) may also transport angular momentum. 

From our numerical results, the dependence of the a parameter associated with local 
turbulence on temperature in a dwarf nova disk can be summarized as follows. The contribution 
of the MRI to the a parameter is significant during outburst. When the temperature is higher 
than the critical value T cr jt ~ 2000 K, the accretion stress is dominated by the Maxwell stress, 
with a ~ 0.1 and 0.01 for the cases with and without net vertical flux, respectively. At quiescence, 
on the other hand, a is very sensitive to the temperature of the disk. If T <J T crit , the stress is 
more than 2 orders of magnitude smaller than that during outburst. 

There are several reasons why global effects may influence these results. Firstly, the presence 
of a net vertical flux can affect the amplitude of the turbulence stress during outburst, that is 
ahot- I n agree with the previous works using the local shearing box (e.g., Hawley et al. 1995; 
1996), our simulations have shown that an sufficient amount of angular momentum transport 
in the hot state requires the existence of a net vertical field. On the other hand, recent global 
simulations obtained a higher stress (a ~ 0.1) even without a net magnetic flux (Stone & Pringle 
2001; Hawley, Balbus, & Stone 2001). There is no observational constraint yet in the structure 
and strength of the magnetic field in dwarf nova disks. Thus, theoretical modeling of the global 
structure of the magnetic field is quite important for a more detailed analysis of a at the hot state. 
Secondly, the ionization fraction of the disk gas, and therefore the importance of non-ideal MHD 
effects, is extremely sensitive to temperature. Thus global models in which the radial and vertical 
temperature profiles in the disk are calculated self-consistently with the local heating rate due to 
turbulent stresses and cooling rate due to radiation are required. Furthermore, at the surface of 
the disk, the radiation from a hot white dwarf could be important source of heating and ionization 
(Hameury, Lasota, & Dubus 1999; Menou 2002). Finally, even when Ohmic dissipation completely 
suppresses the MRI, small amplitude motions remain in the disk and can give a non-zero Reynolds 
stress. In the local model, these motions are damped on a viscous timescale, which may be very 
long. However, they may be affected by global density waves in the disk or tidally induced spiral 
shocks. Thus, global MHD simulations are an important next step for understanding the dynamics 
of dwarf nova disks. 

We would like to thank John Cannizzo, Kristen Menou, and Neal Turner for helpful 
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Observatory of Japan, and the VPP700 at the Subaru Telescope, NAOJ. This work was supported 
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Table 1. Zero Net Flux B z Simulations 
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Table 2. Uniform B z Simulations 
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T[K] 

Fig. 1. — Temperature dependence of the electron abundance n e /n n (dot-dashed curve), the 
magnetic Reynolds number ReM (solid curve), and the Hall parameter X (dashed curve). Typical 
quantities for dwarf nova disks are assumed in the calculations of these curves; n n = 10 18 cm -3 , 
M = Mq, r = 10 10 cm, and va = 10 5 cm s -1 . Arrows denotes the loci of our 4 models for quiescent 
disk (see Tables [l] and |2|) . 




Fig. 2. — Time evolution of the Maxwell stress for the ideal MHD model (dashed curve) and the 
models including ohmic dissipation and the Hall effect (solid curves). Nonideal MHD terms are 
included after 50 orbits, and the size of the terms for each model is indicated in the figure. 
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Fig. 3. — Time- and Volume-averaged Maxwell stress and Reynolds stress as a function of the 
magnetic Reynolds number at 50 orbits {ReM,m) f° r zer ° ne t A ux B z runs. The Hall parameter at 
50 orbits (X50) is also shown at the top of the figure. The time average is taken from 75 to 100 
orbits. 
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Fig. 4. — Time- and Volume-averaged Maxwell stress and Reynolds stress as a function of the 
magnetic Reynolds number at 50 orbits (ReM,5o) f° r uniform B z runs. The Hall parameter at 50 
orbits (X50) is also shown at the top of the figure. The time average is taken from 75 to 100 orbits. 
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Fig. 5. — The critical temperature T cr i t [K], at which the magnetic Reynolds number is unity, is 
shown as a function of the distance from the central star r [cm] for the cases of the neutral density 



n, 



= 10 cm J (circles), 10 (triangles), and 10 (squares). The Alfven speed and the mass of 



the central star are assumed to be va = 10 5 cm s 1 and M = M Q , respectively. 
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Fig. 6. — Time evolution of the Maxwell stress and Reynolds stress for model ZC (ReM,50 = 1 
and X50 = 1). After 100 orbits, the ideal MHD approximation is resumed to mimic a restarted 
outburst. 



